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Abstract 



In problems where the temporal evolution of a nonlinear system cannot be followed, 

a method for studying the fluctuations of spatial patterns has been developed. That 

^ I method is applied to well-known problems in deterministic chaos (the logistic map and 

^v| ■ the Lorenz model) to check its effectiveness in characterizing the dynamical behaviors. 

Q . It is found that the indices fiq are as useful as the Lyapunov exponents in providing a 

CN I quantitative measure of chaos. 

1 Introduction 



An important feature of classical nonlinear systems is that a trajectory traced out by time 
^ I evolution is well defined, so the distance between nearby trajectories is a meaningful function 

^ I of time. The Lyapunov exponents that characterize the distance function have therefore been 

used widely to describe the chaotic behaviors of such systems. Certain quantum systems, 
^ ■ however, do not have such a feature. In particular, self-coupled quantum fields such as 

^ . those in the 0^ theory do not have evolutionary histories that can readily be described by 

trajectories, since the number of degrees of freedom changes with time. In such problems 
alternative criteria for chaos beside the use of Lyapunov exponents must be found. A measure 
useful in the study of QCD parton showers is a set of indices fig that characterize the nature 
of fluctuations of spatial patterns |1[. It is the purpose of this paper to apply that measure 
to classical nonlinear systems and investigate its usefulness as an alternative criterion for 
chaos. 

In microscopic quantum systems it is often impossible to track the time evolution of their 
states without disturbing the systems. Instead, it is the final state that can be measured 
with high accuracy. A prime example of problems of that type is the inelastic collision of 
elementary particles at very high energy, where many particles are created. The momenta 
of all charged particles in the final state can be determined precisely in experiments. Thus 
for each collisional event the momenta of those particles constitute a spatial pattern in 
momentum space. From event to event those patterns change not only in the magnitudes 
and directions of the momentum vectors, but also in the total number of such vectors. 



The challenge has been in finding an efficient way of characterizing the ffuctuation of those 
patterns in experiments where millions of events are measured. Moreover, it has been of 
interest to find out whether the notion of chaos has any meaning for such multiparticle 
production processes. 

In order to answer the latter question, i. e., the meaning of chaos for self-reproducing 
nonclassical systems, it is necessary to apply a chosen measure of fiuctuations in such systems 
to some classical problems for which the criteria for chaos are well known. The issue becomes 
the following: if a classical chaotic system exhibits certain familiar characteristics in its time 
evolution, what can be said about the nature of the spatial patterns associated with its 
trajectories? In finding an answer to this question we shall have accomplished in making 
two beginnings: on the one hand, we shall gain some insight on whether the concept of chaos 
can be generalized to include self-reproducing quantum systems, and on the other, a new 
approach to the study of classical chaos will be opened up. The latter is an unexpected bonus 
that results from the attempts to deal with the demands and concerns of a very different 
field of physics. 

In order to render this paper self-contained, a review of the measure of fiuctuations will 
be given (in Sec. 2) without the particle physics in which it is originated. The body of this 
paper is the application of that measure to the logistic map and the Lorenz attractor [||. 
We compare the dependences of the Lyapunov exponents A on the control parameter r with 
those of the indices /ig. It is the close correspondence between the two measures for both 
deterministic systems that supports our view on the usefulness of /x^ > as a criterion for 
chaos. 

2 Entropy Indices /ig 

Consider the problem of describing a system by making many experimental measurements, 
each of which is called an event. An event consists of a spatial pattern in (i-dimensional 
space. Let Fg be a measure of that pattern to be described below. From event to event Fg 
can fiuctuate. After Af events, a large number, one has a distribution of Fg, which we denote 
by P{Fg), normalized to 1. By taking the normalized moments of P{Fg), defined by 

Cp,, = (F^) / {Fgf , (1) 

we have a quantification of the fiuctuations of the spatial patterns. 

Returning to the definition of Fg itself, it is necessary to recognize first that any descrip- 
tion of a spatial pattern depends on the resolution used. Let the d-dimensional space (call it 
phase space, although it can be just the coordinate space, or just the momentum space, or 
both) be divided into M bins, each having a volume Vbin = S'^- Furthermore, let the intensity 
of the pattern be discretized at the bin level so that at the ^th bin the bin multiplicity 



ni= p (f) rfV (2) 

J Vj 

is rounded out to an integer, where p (r) is the density at the point r. For each event Fq is 



defined, for any integer g > 2, by 

^-^ = ]^ E ^^ (^^ - 1) ■ ■ ■ (^^ - 9 + 1) / ( ^ E riA (3) 

If Qn denotes the distribution of bin multiplicity n in the M bins, normalized to X^n Qn = 1, 



F, = (nM)./« (4) 



then Fq can also be written as 



where n^'^l = n\/{n — q)\ and (■ ■ ■)^ is a (horizontal) average over Qn- By horizontal, we 
mean averaging over the multiplicity distribution in a given event, to be distinguished from 
vertical averaging, such as in (|1|), which is an average over all events. 

The virtue of the normalized factorial moments Fq is that they are trivial for statistical 
fluctuations |0]. Let Q„ be a Poisson transform 

Qn= / -e-''D{s)ds , (5) 

Jo n\ 

where D{s) may be regarded as some dynamical distribution, whose convolution with the 
Poisson distribution (of statistical origin) gives rise to the observed Qn- It is clear that, since 

nM\ = r sW(s)ds , (6) 

' h Jo 

trivial dynamics represented by D{s) = 6{s — n) results m Fq = 1 for all g. Indeed, (^) 
indicates that the statistical fluctuation is filtered out from the factorial moments, yielding 
just the simple moments of the dynamical -D(s). Thus Fq provides an effective description 
of the dynamical fluctuations that generate the spatial pattern under study. 

Now let us consider the nature of the fluctuations from event to event. First, (|l|) can be 
rewritten in the form 

While much information can be revealed by studying all moments p of P{Fq), it is sufficient 
for our purpose here to examine only the neighborhood of p = 1. It is analogous to studying 
the information dimension Di, which is the fractal dimension at order 1 Q. With the 
definition 

E - — r 
" dp "'" 



p=i 



we have, on the one hand, 



E, = (<|.,ln<|.,) . (9) 

On the other hand, if Cp^q has a power-law behavior in M, 

Cp,q oc M'^^(P) , (10) 



which has been referred to as erraticity [Q, then we also have 

InM . (11) 



d 
dp 



Sg OC -i-V'g(p) 



p=l 



For brevity we define 



N = j-^.iP) 



(12) 
p=i 



and refer to them as entropy indices. It then follows that 

in the scaling region, i. e., where Eg exhibits a linear dependence on InM. It is not difficult 
to show how Hq is related to an entropy defined in the event space |l|, || , but that connection 
is not needed here. 

If there is no strict scaling behavior in M, then (p!OD may have to be generalized to 
accommodate a possible scaling law in g{M) 

C,,, ex ^(M)'^^(^) , (14) 



where g{M) is some function of M. In such cases Eg and fig are defined as in (H) and (|T 
the only difference being that M is replaced by g{M) in (jTT|) and (|T^). Thus, instead of ([TT|), 
we would have 

S,(M,r)oc/i,(r)ln(7(M) , (15) 

where we have introduced a control parameter r, the dependence on which has been assumed 
implicitly in the foregoing, but will become explicit in the following sections. The factorizable 
form of (|T3p suggests that g{M) may be determined from J]g{M,r) by evaluating it at a 
particular vq so that 

S,(M,r)oc/3,(r)S,(M,ro) (16) 

where 

Pg{r)=f^g{r)/f^g{ro) (17) 

In this way yUg(r) can be determined only up to an overall factor for all r . 

We have described above a procedure by which one can take A/" events of fiuctuating, 
spatial patterns, and by using (^, (|^), (|^) and (|T^ [or (|TB]) and (|T^] determine a set of 
indices Hq, q = 2, 3, • ■ ■, that can efficiently characterize the nature of the fiuctuations. In 
practice, it is not necessary to examine a large number of fig] /i2 and /is should suffice. In 
the following sections we shall use yU2 as a measure to study the properties of the logistic and 
Lorenz problems and compare its behaviors with those of the Lyapunov exponents A. 



3 The Logistic Map 

The simplest and well-understood example of deterministic chaos is the logistic map 0, ||. 
We consider this example to illustrate the use of /i2, since the value of A for it is well known 
and can therefore readily provide a comparison with our result on /X2- 
In the 1-dimensional interval < x < 1, the map is 

Xj_|_i = rxj (1 — Xj) . (18) 

By repeated iteration one generates a sequence T(xo) = {xq, xi,- ■ ■ ,Xj,- ■ ■}, starting from a 
chosen initial point Xq. Every such sequence can be regarded as a trajectory as time evolves, 
where the time is identified with the number of iterations. The distance dj between two 



trajectories T and T is 

grow exponentially for two near 



Xj Xj 



at the jth step. For r > Vc = 3.5699456 ■ ■ ■, but < 4, dj can 
: Dy trajectories with d^ = |xo — Xq| = e infinitesimally small. 
Except for certain narrow intervals between Vc and 4, A is positive, and the system exhibits 
chaotic behavior. 

The first question to face is how such a behavior in time evolution can be treated from 
the point of view of spatial patterns, which is what fiq are designed to describe. Since a 
trajectory in this case is automatically a collection T(xo) of discrete points in x, the answer 
is, of course, obvious. A judicious choice of a subset of T{xo) is a spatial pattern of interest, 
and each event corresponds to a particular initial value Xq. To see what subset is appropriate, 
we show in Fig. 1 a plot of dj vs j for r = 3.99 and for various small values of dg. The 
value of A can be read off from the initial exponential growth, dj = d^e^^, to be A = 0.66, 
very close to the analytical value ln2 at r = 4. A significant aspect of Fig. 1 is that even 
for do = 10~^^ it takes only 40 time steps for dj to reach 0(1), beyond which dj fluctuates 
with no apparent order. At smaller values of r, but above Vc, A would be smaller and it 
takes longer for dj to get beyond the exponential growth phase. Two spatial patterns having 
infinitesimal do would be nearly the same if the corresponding subsets of T{xq) and T'{xq) 
consist of only the points in the growth phase. To exhibit chaotic behavior it is necessary 
that j > X^Hndo^, so our subset S{xo) C T(xo) should consist of points above that value of 
j. Since we want to study the relationship between A and fi2 for all interesting values of r, 
our choice of points for S{xo) is as follows 

S{xo) = {xa,X2A,--- ,XmA}xg , (19) 

where A = 100, and m = 20. Each event of that type has a specific xq, not included in 
S{xq). We generate A/" = 10^ events whose initial Xq are all randomly generated within a 
small interval (Xo,Xo + 10~^) around an arbitrarily chosen value Xq. For the results to be 
shown below, Xq is 0.35436. Thus all Af events correspond to initially nearby trajectories, 
the distances between any two of which diverge after a certain number of steps. 

For each of the A/" events generated according to the prescription described above, we 
divide the unit interval into M bins of 6 size, count the number of points that fall into 
each bin, and calculate Fg{M) for that event by use of (^. Then Sg(M) is determined by 
performing the appropriate vertical averaging in (^. With focus on g = 2, the dependence 
of E2(M) on InM is shown in Fig. 2(a) for a few representative values of r. Evidently, there 
is no linear dependence. We thus use the generalized scaling form expressed in ([1^ and 



consider the plot of (|I6|). That is done in Fig. 2(b), which shows good hnear behavior. The 
value of tq is chosen to be 3.9. The slopes (32(1^) can be determined from the best fits of 
all the points for each r, and give, by (0), the values of fi2{r) apart from a multiplicative 
constant. 

Figure 3 shows the comparison of A and fi2, where the overall normalization of /i2 in the 
figure is adjusted to agree with A at r = 3.8. The error bars on the values of fi2 are due to 
the deviations from strict straightlines in Fig. 2(b). Clearly, yU2(^) agrees very well with A(r) 
throughout the whole range of r, except that when A(r) < 0,fi2{f') can only be zero, since it 
is a nonnegative quantity. 

It is by virtue of Fig. 3 that we infer the effectiveness of using the positivity of /i2 as a 
criterion for chaos. In fact, fig for higher q have the same property, but they are not needed 
for the simple system under consideration. Thus we conclude that the fluctuations of spatial 
patterns can be used to reveal the chaotic behavior through the study of fig as much as one 
can learn from the temporal evolution of nearby trajectories. 

4 Lorenz Attractor 

We now consider another problem to explore the effectiveness of fig in a dissipative dynamical 
system. The prime example of such systems is the Lorenz model, described by the following 
equations: 

X = -a{x-y) , 

y = rx — y — xz , 

z = —bz + xy . (20) 

We fix, as with Lorenz 0, a = 10 and b = 8/3, and vary r as the control parameter. 
We discretize the time variable and solve (|20| ) by repeated iterations starting from some 
arbitrary point away from the fixed points. The critical value r^ of the control parameter, 
above which the trajectory becomes unstable, depends on the size of the time step 6t used. It 
is found that Vc increases slowly when 6t is decreased. For computational efficiency we have 
chosen 6t = 10~^. Figure 4 shows how rapidly the t dependence of the distance function d{t) 
changes, when r is increased infintesimally from below to above re- We determine the value 
of A from straightline fits of the rising portions of log d{t) for every value of r examined. 
However, because log d{t) does not rise linearly with t for r > r^ a range of values of A can 
be extracted from the fits. We shall indicate the result by shaded bands in A(r). 

We use the same technique as described in Sec. 3 to generate a spatial pattern for each 
event. For r > Vc the trajectory is the familiar Lorenz attractor. Since it is in 3D, we select 
70 points spaced 1 time unit apart (i. e. 10^ time steps of St), and then make a projection of 
them to the x-y plane. Figure 5 shows a typical event. A total of 10^ events are generated, 
each of which starts out initially at a random point in a small cube of size 10"^'^ on each side, 
located at the point xq = 0, yo = 1, zq = 0. Since the Lorenz attractor is confined to a finite 
region of space, which, when projected onto the x-y plane, shows the points mainly along 
the diagonal of x ~ y. We have rotated the coordinates by n/A so that the pattern of points 
is mainly along the new x axis shown in Fig. 5 (—30 < x < 30) with a dispersion in the 



expanded new y-axis (—10 < y < 10). This 2D rectangular space is divided into M square 
bins, and the multiphcity n of points in each bin is counted for the computation of Fg in (^) 
for each event. Using the procedure described in Sec. 2 the quantity J22 is determined and 
plotted against log M in Fig. 6(a) for various values of r. Scaling is obtained by plotting 
against 1^2 ('"o)? as in Fig. 6(b), where tq is chosen to be 28. From the slopes of the lines in 
the latter figure the indices /i2(^) are determined apart from an overall factor, which is fixed 
by normalizing /i2(^) = A(r) at r = 22.9. 

Figure 7 shows the results of our calculations of both A(r) and /i2('^)- As mentioned 
earlier, because of the complicated t dependence of d(t), there is a band of values of A for 
each r. We have determined A(r) only for some representative values of r. Given the errors 
involved, the agreement between A(r) and ^2{r) should be regarded as being quite good. 
The most important point is that they both show stepwise increase at Tc- Thus the utility 
of the positivity of /i2(^) as a criterion for chaos is clearly as effective as that of A(r). 

5 Large M Behavior 

In the previous two examples we determine the slopes l32{r) from Figs. 2 and 6 and by use of 



(p!6|); from /32i'i") we obtain /i2(^) apart from an overall constant. What we want to emphasize 
here is that the scaling behaviors are for a range of M that is not asymptotically large, i.e., 
bin size 6 is not infinitesimally small. For generic problems in statistical physics and fractal 
geometry, the extension toward larger values of M is the conventional procedure. However, 
for problems that we consider here such an extension is inappropriate. To explain that is 
the aim of this section. 

In fractal geometry, for example, one can take the mathematical limit of smaller and 
smaller scale. The fractal object can always be examined with finer and finer resolution. 
But in high-energy physics, on the other hand, the number of particles produced in any 
collisions is finite at finite energy. In the limit 6 —>■ the bin multiplicities can only be and 
1, and all Fg = for q > 2. For the logistic and Lorenz problems that we have examined, we 
have taken a finite number of points (20 and 70 respectively) to display the spatial patterns. 
Thus the M -^ oo limit would also be inappropriate. Knowing exactly where all the points 
are in phase space gives too much information and is not what we seek to determine as the 
measure that can inform us about chaotic behavior. 

What can one say about the large M regions above those considered in Figs. 2 and 6, 
but not large enough to render all Fg = 0? We assert that they are of no dynamical interest. 
For g = 2 it is only necessary to examine the M region in which the bins are small enough 
to contain two or less points in each bin, but not more. Let M^ be the number of bins in 
the eth event with multiplicity n. Then for that event we have 

^2 = ^E^.K-l)/(]^) =2MM,7iV^ , (21) 

where N is the total number of points in the event. If A2 denotes the number of events out 
of the total TV events in which M2 7^ 0, but M„ = for n > 3, then we get 

{F2) = ^,T.MI = 2Mr,{M,)/N' , (22) 



where 



(M2) = ^Y.M^2 , (23) 



M2 



and r2 = A/2/A/' is the fraction of events for which F2 7^ 0, but -Fg%2 = 0- From (0), PT| ) and 
(^) we have 



^l = M'Jr^ (M2) , (24) 



so that from (Rl) follows 



^2 = ^E^eln5e-lnr2 , (25) 

J\l2 e 

when Be = M|/ (M2). In the limit of large M when M2 -^ 1 for nearly all events, then 
Be -^ 1, and 

S2 ~ — lnr2 . (26) 

Now, the probability for a bin in such events to have n = 2 is M~^. Since it can be for any 
of the M bins, we have 



r2 ~ M-^ . (27) 



It then follows from (|T3D that 



fi2 = 1 (28) 

The same line of reasoning leads also to the result 

l^, = q-l . (29) 

In our numerical computation we have verified this result in that Figs. 2(a) and 6(a) exhibit 
straightline behavior at large M with unit slope for all values of r. Since only probabilistic 
arguments have been used to derive the result, it is independent of the structure of the 
model. 

Thus in the search for scaling behavior in problems where N is finite one should not 
go to the extreme large M region just before all Fg — > 0, even though J22 exhibits linear 
dependence on InM there. The behavior that is more relevant to the determination of fig 
involves g{M), defined in (0), as the scaling variable, and it is the plots like Figs. 2(b) and 
6(b) that yield the more pertinent straightline behaviors. 

6 Concluding Remarks 

By working with the two examples, the logistic map and the Lorenz attractor, we have 
demonstrated that the index /i2 is as good as A in marking the chaotic regime of the control 



parameters. One may wonder why the comphcated procedure to determine /i2 should be con- 
sidered when the computation of A is significantly simpler. We reiterate that the rationale 
for studying spatial patterns is rooted in the desire to examine chaotic behaviors in systems 
where following the temporal evolution is not possible, or where trajectories are ill defined. 
Such problems are far more complicated than the simple nonlinear systems considered in de- 
terministic chaos. The complexity of the procedure described in Sec. 2 for the determination 
of the entropy indices fig is commensurate with the complexity of the problems. Applying 
such a tool to study the logistic map seems to be an overkill. But it has to be done in order 
to show the significance of the method. It is only when the agreement between /i2 and A is 
established for problems with known behaviors of A can one claim that /i2 > is an effective 
criterion for chaos, whether the system under study is simple or complex. 

Now that we have a method for treating the fluctuations of spatial patterns, there seems 
to be a wide range of problems that could not be studied effectively previously, but are 
now amenable to analysis by this method. They may range from cracks in dry lake beds 
to galactic distribution. When there is only one event, like the astrophysical problem on 
galaxies, one should divide the whole space into many subspaces, each constituting an event, 
study the multiplicity fluctuations in bins of various sizes in each subspace (event) and 
then average the fluctuations of those patterns over all subspaces. Even in problems of 
conventional deterministic chaos, it is not always easy to flne-tune the initial conditions 
experimentally. Studying the properties of spatial patterns may allow an experimentalist 
to circumvent the flne-tuning difficulty. It would be very interesting to explore through the 
study of fi2 the possible universality among many flelds that have hitherto been regarded as 
totally unrelated. 
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Figure Caption 

Fig. 1 Exponential growth of the distance dj between two trajectories as the time step j is 
increased. 

Fig. 2 Behaviors of S2 for the logistic map as a function of (a) InM and (b) S2 {tq) for 
various values of the control parameter r. The value of tq is chosen to be 3.9. 

Fig. 3 A comparison of ^2 with the Lyapunov exponent A for the logistic map. 

Fig. 4 The behaviors of the distance function d{t) for the Lorenz attractor at two values of 
r close to r^. 

Fig. 5 The spatial pattern of one event for the Lorenz attractor when projected onto the 
x-y plane and rotated by 7r/4. 

Fig. 6 Same as for Fig. 2, but for the Lorenz attractor, and with tq = 28. 

Fig. 7 A comparison of ^2 with the Lyapunov exponent A for the Lorenz attractor. 
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